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Firstly, we present a reformulation of the standard canonical approach to spherically symmetric 
systems in which the radial gauge is imposed. This is done via the gauge unfixing technique, which 
serves as the exposition in the context of the radial gauge. Secondly, we apply the same techniques to 
the full theory, without assuming spherical symmetry, resulting in a reduced phase space description 
of general relativity. The canonical structure of the theory is analyzed. 
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I. INTRODUCTION 

General relativity viewed from the canonical perspective is a theory governed by constraints. The 
constraints are first class [l| and, therefore, generate gauge transformations. These transformations, 
being the canonical realization of the general principle of relativity, represent not only the invariance 
of the physical content of the theory with respect to changes of coordinates, but are also intricately 
related to the dynamics of the theory. In fact, general relativity is a fully constrained theory, meaning 
that the Hamiltonian consists only of constraints, which is the central issue in the so-called problem of 
time. One of the ajmroaches introduced in the literature to deal with the constraints is the procedure of 
deparametrization [2,I3|. Usually, the procedure makes use of some carefully tailored matter content of 
the theory to define physical coordinates and, therefore, fixes the freedom in the choice of coordinates. 
A recent application of deparamterization was discussed in Q, where it has been used in two ways. 
Firstly, the presence of nonrotating dust was exploited. Deparametrization with respect to it endowed 
the theory with a raeferred notion of time. Such a construction had been known in the literature Q- 
The new input of was to deparametrize the remaining gauge freedom, of spatial diffeomorphisms, 
with the use of certain geometrical quantities. The construction was based on the introduction of an 
observer: a point in the spatial manifold, and a preferred frame which represented the spatial directions 
as described by the observer. It resulted in a construction of observables invariant with respect to spatial 
diffeomorphisms (preserving the observer). Moreover, the paper contained a derivation of the Poisson 
algebra of the observables and some remarks on the possibility of using that construction for solving the 
diffeomorphism constraint and, hence, obtaining a reduction of the phase space of general relativity. Such 
a reduction is completed in the current paper. 

A related construction of observables has recently been of interest , partially in the context of the 
AdS/CFT correspondence. The main difference with the construction presented in Q and underlying the 
current phase space reduction is that in those works, the Hamiltonian constraint is also being gauge-hxed 
geometrically. By the virtue of the additional gauge condition, the geodesics involved in specifying the 
observables are spacetime geodesics, as opposed to spatial geodesics in Q. 


II. MOTIVATION 

The construction of the observables in Q can be viewed as expressing the canonical fields in a certain, 
carefully chosen, geometrical coordinate system. It is a system of coordinates centered in a point which 
represents the location of the observer. Points of the spatial manifold are then labeled by their proper 
distance from the observer and the angles which describe the direction in which a given point is located 
with respect to the observer. The usage of the proper distance from the observer is a very peculiar feature 
of the construction, and it is transferred to the reduced variables suggested by 0- Namely, it implies 
a special role played by the surfaces of constant radial distance from the observer. It turns out that 
the internal geometry of those surfaces (the components of the metric tensor tangent to them) and the 
tangential components of the momentum of the metric are the variables which suffice (together with fields 
describing the matter content of the theory, in case they are present) for the description of the geometry. 
On one hand, this fact enables the reduction of the description of the system to those variables. On the 
other hand, it seems naturally suited for describing physical systems whose configurations favor a system 
of concentric structures. Such a situation takes place in the case of spherically symmetric configurations. 
There the concentric surfaces become spheres, and their internal geometry can be described by a single 
function related to the area of a given sphere. 

The current paper presents a derivation of the reduced phase space, which stems from the insights of 
j^. The derivation follows a route of gauge fixing, and since spherical coordinates similar to the ones 
used here have been used many times in the literature, we chose to call the choice of the gauge we are 
using the radial gauge (see, e.g., Q). It leads to reduced variables that have Dirac brackets identical with 
the corresponding Poisson brackets of the observables constructed in Q. The derivation concerning the 
full theory is preceded by one performed in the context of spherical symmetry. This serves two purposes. 
Firstly, it exposes various features of the radial gauge and of the methods that are employed in the current 



3 


paper, e.g., in the spherically symmetric setting, a certain physical choice is more transparent and the 
nonlocalities which appear in the Hamiltonian are easier to understand. Secondly, the peculiar feature of 
the reduced variables of the full theory opened a new possibility for a quantization of the theory. After a 
translation to connection type of variables and a quantization relying on loop quantum gravity techniques, 
the geometry becomes effectively 2 + 1 dimensional, which is a significant simplification. Because of the 
peculiar feature of the radial gauge, the quantum theory obtained in this way can readily accommodate 
a quantum definition of spherically symmetric configurations. In fact, there is more then one definition 
possessing different advantageous properties. A detailed discussion of the quantization of the reduced 
phase space obtained in the current paper, together with the identihcation of the spherically symmetric 
sector of the theory is presented in a companion paper and summarized in Ea¬ 


rn. SPHERICALLY SYMMETRIC MIDISUPERSPACE 

In this section, we will apply the radial gauge to the spherically symmetric midisuperspace treatment 
of general relativity. On one hand, this serves as an introductory exercise to better understand this 
gauge in the context of full general relativity and to compare our results to previous work in the context 
of midisuperspaces. On the other hand, we use the resulting Hamiltonian as an ingredient in one of 
the definitions of spherical symmetry in the companion paper @ (see also El), since it preserves the 
respective quantum states as an operator on the full Hilbert space. 


A. Canonical strncture 

We start following the treatment of the Arnowitt-Deser-Misner (ADM) formalism for spherical 
symmetry performed in EH- The coupling of nonrotating dust is discussed in |l3| (and references therein). 
We will briefly recall the main steps here. 

The gravitational part of the action we will work with is^ 

= ^J dtd^aN^^{K^^K,, - . (I) 

The spherically symmetric midisuperspace sector of the ADM formulation of general relativity can be 
obtained by restricting the spatial line element to 

ds^ = K^{r,t)dr'^ + R^{r,t)d^^, ( 2 ) 

while the restriction affects the lapse function N{r,t) and the shift vector held N{r,t) by restricting the 
shift to only having a radial component N''{r,t). Both of them are, of course, only functions of r and t 
in the symmetric context. The canonical analysis performed in the given variables leads to the following 
Poisson brackets 


{R{r), PR{r)} = 6{r-r), {A{r), PA{r)} = S{r - r), (3) 

where the momenta are dehned by 

PR = -^{R-N-R')-^{A-{ANy), P^ = -^{R-N-R'), (4) 

where a prime denotes a radial derivative and a dot, a derivative in time. The spatial diffeomorphism 
constraint retains only its radial component and reads 

C[iV] = drN^PRR'-APj^ + C^^^^), (5) 


^ Notice that to restore the units convention G = 1 = c used in [I2ll . one should put x - ^ restore the units convention 

commonly used in quantum gravity literature, namely, SttG = 1 = c, one should put x = Stt. Note also that the latter 
convention will be used throughout Section HYi 
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whereas the Hamiltonian constraint is given by 


J f" oo 

drNh = 

0 



PrPk 

R 


) + x( 


RR" 

A 


RR'K' R'^ A, 
A2 2K 2’ 


( 6 ) 


where and are contributions from matter fields (in case they are present in the theory). 

The addition of nonrotating dust to the theory introduces terms into both of the above constraints; 
however, after deparametrizing with respect to the dust field (following the prescriptions discussed in 
i) , the vector constraint regains the above form, while the Hamiltonian constraint is replaced by a true 
Hamiltonian equal to 


J f" 00 
0 


drh. 


( 7 ) 


B. Mass functional 


As it was noted in (T^l, in spherical symmetry, it is possible to give a local definition of a mass functional. 
One can check by explicit calculation that the expression 


m := 



( 8 ) 


has weakly vanishing Poisson brackets with the vacuum vector and Hamiltonian constraints. Moreover, 
its spatial derivative vanishes on the constraint surface, and therefore, the functional is in fact a constant 
throughout space and time. The reason why this functional represents the mass becomes apparent once 
asymptotic conditions and possible boundary terms are taken into account (this is discussed in Appendix 
[Bt or when the Schwarzschild solution is investigated (see Appendix [C]) . 


C. Gauge (un)flxing — measuring distance from zero 

The next step is to implement the radial gauge. We choose 

A(r) = 1, (9) 

since it corresponds to choosing the radial coordinate r to measure the proper distance. In this section, 
we will discuss a construction which corresponds to having r measuring the distance from the center of 
symmetry (or zero). In Section IlIID1 we will discuss an alternative construction in which the distance 
will be measured “from infinity”. 

A technique equivalent to employing the Dirac bracket (with respect to the gauge fixing ([H]) and the 
spatial diffeomorphism constraint), is to use gauge unfixing (l3. [l5|. This amounts, in the first step, to 
modify the momentum P\ in all constraints except for C by powers of C in such a way that it Poisson 
commutes with the gauge fixing condition on the gauge fixing surface, thus having equal Dirac and Poisson 
brackets. In a second step, one can interpret C as a gauge fixing for © and drop it. Since © is Poisson 
commuting with the remaining constraints (in this case, only the Hamiltonian constraint) after P\ has 
been modified as described above inside the other constraints, becomes a first class constraint. As 
said before, this is classically equivalent to employing the Dirac bracket or using ([H]) to gauge fix C^. In 
this paper however, we will not drop C from the list of constraints, but merely use techniques developed 
in the context of gauge unfixing [l5j to express the Hamiltonian in terms of the reduced variables, see 
equations (rT5]l and (fTTll (in other words, to compute a gauge invariant extension thereof^). 


2 Changing the constraint structure, however, might be advantageous in the quantum theory, which was the original 
motivation for introducing gauge unfixing in Q- 

^ Practically, this means adding a power series in the vector constraint to P\ and inserting the resulting expression in the 
Hamiltonian. In the case of general relativity (coupled to matter fields), the series for Pa terminates after the first term 
and one obtains a manageable expression for the Hamiltonian. 
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There are two main reasons for which we choose to use the gauge unfixing method. Firstly, it can be 
seen as a method of directly finding the Hamiltonian preserving the gauge, which in some cases (e.g., 
the nonsymmetric case in this manuscript) leads to a more straightforward way of obtaining such a 
Hamiltonian. Secondly, it provides a clear phase space picture of the process of implementing the gauge 
fixing (see Fig. [T]). However, as has been said before, the gauge unfixing method is equivalent to the 
standard Dirac procedure of implementing the gauge fixing. This can be best seen in sections IlII El and 
IIV D1 where the standard approach of solving the constraints in a given gauge is discussed. Also, the 
shift vectors discussed in Appendix |A] are closely related to the shift vectors needed to invert the Dirac 
matrix if one chooses to use the Dirac bracket approach. 



FIG. 1: The picture illustrates schematically the gauge unfixing procedure. F is a phase space in which 
the theory is defined. The lined, red (for colors see online) surface C is the constraint surface on which 
the lines represent gauge orbits. The vertical, blue surface G is the surface on which the gauge fixing 
condition is satisfied. The thick black line D = C nG denotes the points in F that we are most 
interested in. The procedure of gauge unfixing provides us with an extension of the Hamiltonian off C 
such that the Hamiltonian vector field of the extended Hamiltonian is tangent to both C and G, and 
therefore also to D. In fact, gauge unfixing defines dynamics (equivalent to the one taking place in G) 
on all of G; however, it is the restriction to D which allows for a simple link between the two. 

Consider an expression linear in P/^, that is, the field smeared with a phase-space-independent function 
p. Our goal is to find a function M(f) such that, denoting 

J f" OO oo 

drPA{r)fi{r) = / drPA{r)fi{r) - C[M] (10) 

0 Jo M 

the tilded expression satisfies 

{A(r), f drPAir)^i{r)} = 0 (11) 

Jo 

on the gauge fixing surface A(r) = 1. This condition yields the equation 

drM{r) = ( 12 ) 

iJ- 

Since M, being a smearing field of the vector constraint in the spherically symmetric context, has to 
vanish at zero,^ we have 

M(f) = f df^{f). (13) 

M Jo 


^ This is the case for regular spacetimes. See Appendix |C] for a discussion of an example of a more general setting. 
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It will be convenient to have an expression for P\(r)^ which we can obtain from the above procedure, 
choosing /x(f) = 6{f - f). It follows that 


PA{r) = PA{r) - J drCr{r). (14) 

Notice that on the gauge fixing surface 

PA{f) = PA{f) - r df{PRR' - P; + (r) (15a) 

= PA(f)- Pdf(PjiP' + C“"“)(f)-PA(f) (15b) 

Jr 

X oo 

df(PflP' + C“"“)(f), (15c) 


where we used the fact that Pa(?') ^ 0 when r oo.® It is not surprising that Pa drops out from PA(r), 
since by the virtue of m on the gauge fixing slice 


SPa 

SPa 


= 0 . 


(16) 


We are now ready to gauge (un)fix our model. This is done by setting A(r) = 1 and using the expression 
(I15c|) in the Hamiltonian. The Hamiltonian turns out to be 


^[^]lgauge-fix 

■ X" (X"^ X~‘'d^»«'" t"”") (") 


+ + (17) 


It is instructive to look at the procedure of gauge fixing applied above from the perspective of the 
unreduced phase space. As it is argued in Appendix lA 11 employing the radial gauge amounts to fixing 
the shift vector to be 

^Hir) = ^f^dfN{f)^^{f)-^{f) + ^{f)J_ drCrir)y (18) 

In case one couples the model to nonrotating dust and deparametrizes the Hamiltonian constraint with 
respect to the proper dust time, one obtains the true Hamiltonian in the form 


^truel 

I gauge-fix 

■ X" * ( 2 XiJ (X" <">)'" p X" - 5)) ■ 

To sum up, we have solved the spatial diffeomorphism constraint via the gauge unfixing procedure, 
using the gauge A(r) = 1. The reduced phase space turns out to be coordinatized by R{r) and PR(r), 
subject to the Hamiltonian constraint (HB. Since the Hamiltonian constraint Poisson commutes with the 
gauge fixing condition due to the addition of spatial diffeomorphism constraints, the Dirac and Poisson 
brackets are equivalent, and the dynamics can be readily computed using either one. 

Having obtained the reduced phase space, let us compute the equations of motion for the variables in it. 
For the sake of conciseness, let us spell them out for the case deparametrized with the use of nonrotating 


^ This condition is satisfied for spacetimes which are asymptotically flat. See Appendix HI for a more detailed discussion of 
asymptotical flatness. Note, that in the case of the construction from Appendix o this condition is not satisfied; however, 
in the construction presented there, it is also not necessary. 
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dust (with no additional matter) and using the convention x = 1- They turn out to be 


R{r) 


Mr) 





Mr) 

R{r) 


^ PR{r)F{r) 

R^{r) 


Fjr) \ 

R^{f)) 
F^{r) 
R^{r) 


+ Pkir) 



Fjr) \ 

i?2(f) ) ’ 


(20a) 

(20b) 


where F{r) = - df {PrR') (f) is the expression (jl5c|l of Pa in terms of the reduced variables. These 
equations are hard to solve in general; however, we can readily identify the trivial solution. Putting 
Pnir) = 0 constantly in time we find 

R(r) = 0, (21a) 

0 = i?"(r), (21b) 

which has a unique solution R{r) = r (since by construction i?(0) = 0 and R'{0) = 1), which is a relation 
characteristic for flat spacetime. Note that the Hamiltonian vanishes for this solution. As has been 
demonstrated in sm, Il7l |. the Hamiltonian is equal to the momentum of the nonrotating dust field we 
used for deparametrization. That momentum, in turn, is proportional to the rest mass density of the 
dust. Therefore, the vanishing of all of those quantities represents a limit in which the dust is not present 
in the theory, and the spacetime is just the Minkowski spacetime. 


D. Gauge (un)flxing — measuring distance from infinity 

In this section, we will present a gauge fixing procedure closely related to the one presented in Section 
IIII Cl The difference here is that we want our gauge fixing to correspond to setting r measuring the proper 
distance “from infinity”. 

The gauge fixing condition we use in this section is 

A(r) = 1, (22) 

like it was previously. Finding a gauge invariant extension of Pa (see equation (1101) 1 amounts to solving 
the equation (see equation (1121) 1 

drM(r) = fi(r). (23) 

Instead of integrating it from zero, as we did in the previous section, here we will use the fact that the 
shift vectors we consider vanish at infinity, which leads to 

M{f) = - / dffi{f). (24) 

M Jr 

This formula leads to the following expression for Pa as a function of the reduced phase space variables 
(compare with equation (I15cll l 

pA{f) = £ df{PRR^ + C^^^^){f), (25) 

where vanishing of Pa at zero has been used. The Hamiltonian of the reduced phase space is (compare 
with equation (ITO)) ! 

-^[-^llgauge-fix 

■ X‘ g . C"-) (f) 

+ x(RR"+^-i)+(.“'j. ( 26 ) 

It corresponds to fixing the the shift vector to be (derivation is presented in Appendix lA 21) 


( 27 ) 













E. Where should we measure from? 


In this section, we will compare the two alternative variants of implementing the radial gauge, discussing 
the differences and highlighting the advantages. 

The difference between the two ways of “measuring” lies in the domains over which certain quantities 
are integrated. We use the terminology “measuring from zero” because the expression for T’a(?') spelled 
out in equation (I15cl) generates transformations of the fields in the reduced phase space in the domain 
from r to oo. This is what is expected if r measures the distance from zero, since then a variation of 
the metric at a given r should affect the phase space variables in the radial gauge in points lying further 
than r. This fact becomes straightforward when translated to the language of the observables defined in 
[3|, since the variation of the metric at a given point affects the radial geodesics used in the definition 
of the coordinates needed for the observables only after they cross the perturbation. Hence, it is the 
“measuring from zero” choice which corresponds more closely to the construction of observables from Q. 
On the other hand, “measuring from infinity” might be more natural in some settings. Firstly, as we 
will see in Section IIV Ci the option which was suggested in the last section of [3| for a solution of the 
vector constraint corresponds to Secondly, it naturally accommodates the Schwarzschild solution 

as discussed in Appendix 0 

It should be noted that the two expressions for Pa are closely related. The vector constraint in the 
radial gauge can be rewritten in the form 

P;(r) = (Pfli?' + C'”"“)(r). (28) 

Integrating it from 0 to oo, using the fact that Pa(0) = 0 and assuming® limj.^oc. PA{r) = 0, the following 
condition is obtained 

+ = 0. (29) 

This means, that if asymtoptical flatness (in the sense of Appendix [B]) is assumed, this condition needs 
to be satisfied throughout the evolution of the system. However, in such a case, the left-hand side of 
Poisson commutes with the Hamiltonian, and therefore, it can be viewed as a restriction on the initial 
data consistent with the dynamics. 

The condition (1291) renders the two expressions for Pa numerically equivalent (they necessarily sum up 
to zero). However, they are distinct when viewed as functionals on the reduced phase space of the theory, 
since they generate diffeomorphisms of the canonical data in different domains. 


F. Summary of the spherically symmetric setting 

In the context of canonical spherically symmetric general relativity, we have presented derivations of two 
reduced phase spaces. Each is parametrized by two functions, R{r) and Pii{r), of which the first describes 
the area (given by 47ri?^) of a sphere located at the proper distance r from the center of symmetry, while 
the latter is its conjugate momentum. Their dynamics is governed by the Hamiltonian constraint (ED or 
(I26p (in case we additionally deparametrize with the use of nonrotating dust, by the true Hamiltonian 
(fIT?D or its analogue), yielding equations of motion spelled out in ((^01) or their analogue. The reduced 
phase spaces support asymptotically flat spacetimes (see section IB II of Appendix |B] for details), and the 
initial data is subject to the condition (1^91) . One of the two phase spaces may be modified such that the 
Schwarzschild solution can be readily identified, for details, see Appendix [C] 


It is a consequence of the asymptotical flatness conditions discussed in Appendix m It is is not true, however, in the 
generalized setting of Appendix [Cl 
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IV. FULL THEORY 
A. Canonical strncture and gauge fixing 

We will now turn to the full theory and rederive the main results of [il via gauge fixing. This way 
of derivation is complementary to the explicit construction of the observables as performed in Q and 
leads, as is clear by general arguments [^, to the same observable algebra. It was the basic idea of 
to define a physical coordinate system by specifying a point (Tq in a spatial slice S of the spacetime 
manifold and using the exponential map from the tangent space to E to define such coordinates. 

Since the exponential map relies on radial geodesics, such a choice of coordinates can be referred to as the 
radial gauge. The basic concept of using such a radial gauge is, of course, not new; however, the detailed 
description of the observables and their Poisson algebra has only been worked out recently in . 

In order to define local (Gaufi) coordinates (x^), 1=1,2,3, we need to choose three linearly independent 
vectors in To-pE which transform properly under spatial diffeomorphisms. This can be accomplished as 
follows: at ctoj fix a frame eo/, where, in general, eo\eQj5^^ + and f, j,... are spatial tensor indices on 
E in some general coordinate system. Given we now construct the frame e“ with as 

e\=Y, Mijeo'j, (30) 

J=i 

where Mjj is a lower triangular matrix.^ Now, given an infinitesimal diffeomorphism acting on , we 
have Qij ->• Qij +£]\[qij. The frame e\ constructed as above could a priori transform as e\ CNe\ + Lj'’ 
for some antisymmetric matrix Ljj. However, it follows from (1301) that Ljj = 0, meaning that the frame 
e\ transforms exclusively under the same spatial diffeomorphism that g^ transforms under. From this, 
it follows that the map exp(x^e)') is spatially diffeomorphism invariant for spatial diffeomorphisms 
vk satisfying 


4'(cro) = cro, '^'{ao) = M (31) 

for some lower triangular matrix M. In coordinates, the second condition translates to djN''= 
e}di{N^ej){ao) on the vector field iV* generating ^k. Such diffeomorphisms will be later referred to by 
Diffobs as they where coined in Q- These coordinates have only a finite range in general; however, they 
are always defined in a neighborhood of uoi and their range is maximal in case of spherical symmetry (up 
to nontrivial topology). In what follows, we assume they are defined globally. 

Given the Cartesian “Gaufi” coordinates (x^), we can introduce spherical coordinates (y“) = {r,9). r 
then labels the geodesic distance from ao, and 9 collectively denotes the two coordinates on surfaces of 
constant r. We will denote the two angular coordinates contained in 9 as A, B ,and a = r, A, following 
the notation of (J. Tensor indices can now also be written in terms of these coordinates. For a general 
tensor density T of weight w and general coordinates ( 2 *), we can write both 


rjil J... 

^ KL. 


cd.. 


d{x^,x^,x^) 

d{z\z^,z^) 

5(r,6»i,6»2) 


dx^ dx’^ dz^ dz^ j,- 

dyO. Qyb Q^k Q^l 

dz'- dz^ dy'^ dy'^ ' 


(32a) 

(32b) 


We will call the two introduced coordinate systems adapted to the metric qtj. Note that although both 
(x^) and (y“) depend on the metric qij, the relation between them {x^ = rn^, where is a unit vector 
field depending only on 9) is independent of the metric. 

We will work within the ADM formulation [a of general relativity coupled to matter fields, that is 
with canonical brackets 


{<7y(cr), /'(o-)} = (5(cr,cr)(5(V'), 


7r“ (cr)} = S{a,a)S2 ■ 


( 33 ) 


In fact, Cj is the result of the Gram-Schmidt orthonormalization process performed on eoj with respect to qij. 


7 
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i,j,... are spatial tensor indices on S and a, a' label schematically the different matter field species. We 
restrict here to scalar fields for concreteness. 

Basically, one would now like to introduce the gauge q^a = where r,a are coordinates adapted 
to Qab, in order to use the coordinates a = r, A as the physical spherical coordinates described earlier. 
However, this gauge condition does not impose any restriction on since any metric satisfies this 
property when expressed in its adapted coordinates. We therefore specify a reference metric qij on E. 
This metric induces the spherical a and the Cartesian I adapted coordinates via the above procedure. 
By construction, we have qfa{<j) = Sfa- It will now be crucial to observe that given a second metric qij on 
E, which in the coordinates adapted to qij reads qfa = Sfa, we have qfj,dy°‘di/' = qfj,dy°'d'l/’ + q^gdy^dy^ 
for some q^jjj- In simple terms, q^i and q^^ differ only in the AB components, if and only if, they induce 
the same adapted coordinate system Q. Moreover, if qfa = Sfa, then the coordinates adapted to qij and 
qij coincide. 

Using the coordinates adapted to qij, we can impose the gauge condition qfa{<j) = Sfa for the spatial 
diffeomorphism constraint. Since the “check” and “hat” were introduced only for clarity of the argument 
in the previous paragraph, we will now stick to the “checked” coordinates dropping the notation with the 
“check”. To proceed with the gauge fixing, we need to show that the Poisson bracket of the gauge fixing 
condition with the vector constraint, namely, 

{g™(a), C[iV]} = 27V(,.,)(a) (34) 

is invertible on the constraint surface. This problem has already been solved in Q, where the equation 

2V(,.A^a) = Wra (35) 

for arbitrary (symmetric) ojra has been solved for the vector field generating an element of Diffobs- 
The solution reads 


N{r,e) = 


^-ujKj{Q)h'^^rn^ 




/ dr U)rr{'f'^0) 

Jo 

df q^^{f,9)^UJrAif,d) dr UJrr{r,0)^^ 


dB, (36) 


where is a unit vector field such that = rn^, = S^'’ - n^n'^ and lukj is built from the elements 

of liJkj in the following way® 


"Oil 

012 

Wl 3 


LOii 

0 

0 ■ 

021 

UJ22 

W23 


OJ2I + W12 

UJ22 

0 

031 

UJ32 

W33. 


tcai + wia 

UJ32 + W23 

W33. 


In the following, it will be convenient to rewrite this solution as® 

= f^d^dD-^‘^\a,a)ujrb{d) 


with D ^“^(cr, cr) satisfying 


f d^aDab{a-,a)D = SaS{a,a). 

JY. 


(37) 


(38) 


(39) 


One more useful notation we introduce is the smeared gauge fixing condition, namely, for a field v°' 


Qra 


[i/“] := J^d^aqra{<j)i^°‘{cr)- 


(40) 


® The bar in uJ, being a notation inherited from 0, denotes a certain operation on uj, unlike in f and a, where it is used 
to denote a point different than r or a, respectively. 

® Note, that the integrals in equations from ll38l l to gg are not really geometrically defined, as the notation may suggest, 
fn fact they are defined as integrals in the “checked” coordinates mentioned above. 
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Note that with the notation we introduced, we can rewrite (IMl) in the form 

C[N]} = f^d^a f^d^aiy‘^ia)Dabi<7,a)N\a) =: D[^,N]. (41) 


The Dirac matrix reads 


{ 


( "1 ( C\m a n -D[u,M] \ 


(42) 


and can be easily inverted on the constraint surface (^[-/Vobs] = 0 = q^a - dra by using (I5IJ1) . The Dirac 
bracket between two phase space functions F and G now reads 


{F, G}^^ = {F, G}- J^d^ a J^d^a{F, C,{a)}D-^<^\a,a){q,t{d), G} 
- f^d^a f^d^a{F, qra{<j)}D-^‘^\a,a){Cbid), G}. 


(43) 


Let us now compute the Dirac brackets between the elementary phase space coordinates qAB,P'^^, Qra,P^°' 
to test if we recover the same algebra as the corresponding observables formed in Q (see also remarks 
on that algebra discussed in HI)- First, the bracket 

{qABicr), p^^id)}^^ = {qABicr), P^^id)} = 6{a,d)SfAdB) (44) 

follows directly. Next, we compute 

{qra[v°^], ^’Idb = F}-J^d^a J^d^d{qra[i^°-], Cc{cr)}D~^‘'’’{a,d){qrb{d), F} 

= {qra[’^°'], F} - J^d^a J^d^d J^d^diy°'{d)Dac{^,(7)D^^‘^^{a,d){qrbid), F} 

= {qra[iy‘^l F} - {qra[l^^l F} = 0 (45) 


for an arbitrary F, which exemplifies that it is consistent to impose qra = 5ra before computing the Dirac 
bracket. This corresponds to the identity Qraii’-,d) = dra for the diffeomorphism invariant observables of 
j^. The most interesting Poisson brackets involve p’’“. For F independent of p^^, we compute (for clarity, 
we consider p^°' smeared with a field Ka) 

{F, p™M}db = {^. p’^^lKa]}-fj^a£d^d{F, Gb(a)}D-^^^(a,d){qrc(d), 

= {F, f^d^a f^d^d{F, Gb(a)}D-^’’^(a,d)K,(d) 

= {F, (46) 

Also, this Dirac bracket agrees with the observable algebra from Q. Finally, the bracket 
{p [Ka], follows analogously. 


B. Gauge (un)fixing in the full theory — measuring distance from zero 

Complicated Dirac brackets such as (HSI) are an obstacle to quantization. While the reduced phase space 
is already coordinatized by the canonical pair qAB, P'^^, as well as the matter fields, the Hamiltonian 
still contains p'^°‘ as a problematic part. Therefore, one either needs a representation of the Dirac algebra 
containing p’’“, or we need to express p’'“ classically in terms of the remaining variables by adding 
constraints. A means to do this is the formalism of gauge unfixing [alii implemented above in the 
spherically symmetric case. Here, it boils down to computing a gauge invariant extension of a phase 
space function canonically conjugate to the variable being gauge fixed, namely p'’“, with respect to the 
gauge flow of the gauge fixing condition “ dr a = 0, by adding terms proportional to the Diffobs spatial 
diffeomorphism constraint. Let us consider a phase space independent, symmetric smearing tensor pab, 
which will later be limited to have only some nonvanishing components. We are interested in 

J P°-^Pab = J P°'^Pab- G[M], 


( 47 ) 
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where the vector constraint C contains both the gravitational part and possibly a contribution from some 
matter content of the theory (denoted below by The vector field M should be chosen such that 

A* 

the condition 

[qraj J = 0 (48) 

holds on the gauge fixing surface. This condition translates to 

Mra = 2M(r;a), (49) 

fj. 

which is an equation of the type (1351) and therefore, we can readily spell out its solution 


M(r, 0) 


Jq dr flrr{r,9) dr 

df q^^{f,0)^fj.rA{r,9) 




Ob- 


(50) 


For the time being, let us consider only fields fi which are vanishing at zero. For such fields, the first 
term in the above solution does not contribute. Considering first /r to be such that 


PrA{r,0) = S{r,ro)S{9,9o)qAAoir,d), (51) 

we obtain 

M(r,9) = Q{r-ro)S{9,0o)dAo- (52) 

When implemented in ([Tf)) , the above vector helds give 

lfAoiro,9o) = £ dr^VBP^Aoir,9o) - 0o)j + Jim/AoC?-,^* 0 ), (53) 

where a notation from (J has been employed, namely, D denotes the derivative covariant with respect to 
the tangential metric qAB, and the zero above the symbol means the derivative is with respect to the 9q 
variable. Until the end of Section llVBi we will drop the last term from the above formula, assuming it 
is zero. As will be discussed later (see Section PV PI and Appendix IB 21) . this is problematic.^'' However, 
for the sake of clarity of exposition, we adopt this assumption for the time being. Note that under that 
assumption, the expression on the right-hand side does not depend on which is consistent with the 
condition (1^ . 

Secondly, let us consider ^ to be such that only the 


firrir,9) = 6{r,ro)S{9,9o) 


(54) 


component is nonvanishing. It gives 

M{r,9) = ^Q{r-ro)S{9,9o)dr-^Q{r-ro) f dfq^^{f,9){dAS{9,9o))dB, (55) 

M Z Z Jro 

which leads to the following form of (1471) 

1 ^00 

f r{ro,9o) = -- / dr((p^^gAiJ,r-)(r,6»o) + C'“®'**(r,6»o)) 

Z Jro 

X oo 0 / ^00 / 0 1 W 

^ drVAh^^M) £ df\VcP^B{f,9o)--Cr'\f,0o)jl (56) 


In fact, to be mathematically precise, the limit on the right-hand side of equation 15^ should be written in front of both 
terms together, because only together, they have a finite limit. 
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where an assumption analogues to the one simplifying equation (1531) was adopted. Note, that expressions 
(I53p and (I56p realize our aim; namely, they provide expressions for radial components of the momentum 
p°‘^ in terms of the reduced variables: the purely angular components of the metric and its momentum 
and possibly the matter fields. 

The only part of a we have not addressed yet is its behavior at zero. However, since is a tensor 
density, the components p^a vanish at zero. Therefore, we choose jfa to also vanish at zero (since it has 
to do so on the constraint surface anyway). One might worry that in this way, we disregarded the first 
term in ([50]), but we will see later that such a choice is dynamically consistent. Actually, to guarantee 
that indeed jP~a given by (1531) and (1561) vanish at zero, we need to impose conditions on the canonical 
data analogous to the condition (l29l) . namely 

r°° / 0 1 \ 

driVBP^Ao{r,Oo)--CZ'\r,0o)\ = O, (57a) 

1 r 

~2 L + C^r“(A0o)) 

r 0 / r oo / 0 1 W 

+ drVAh^^{r,eo) df\VcP^B{f,9o)0o)\Uo. (57b) 

Fortunately, also in the current case, these conditions turn out to be preserved by the dynamics, and 
hence, are just restrictions of the initial data. 

Knowing p^a, we are now ready to spell out the Hamiltonian for the reduced phase space. It is 

^[^]lgauge-fix = / drd^ 0 Ni^-^G - + /,-tt j ^ (53) 

where 

G = ^{p^r)^ + 2q^^p^ aVB - qABP^^Vr + {qACQBD - \qABqcD)p^^P^°, 

• q^^qAB,rT - \q^^,rqAB,r - \{q'^^qAB,TY, ( 59 ) 

detg = qeeq<j„j, - (qe^^- 

The shift vector that corresponds to this Hamiltonian is presented in Appendix IA3I 


C. Gauge (un)fixing in the full theory — measuring distance from infinity 

Also in the full theory case, an alternative implementation of the radial gauge can be considered. The 
vector constraint in the radial gauge has the form 

drP\ = -Vbp^a + \cT\ (60a) 

drP\ = \qAB,rP^^ " 5^. + (60b) 

The expressions for the ff a present in ( 1551 ) and ( 1551 ) correspond to integrating the above equations from 
infinity. As it was pointed out in j3|, the components p^ a of the momentum have to vanish at zero due 
to its tensor density character. Therefore, we can integrate the above equations to obtain (we use the 
tilde to underline that it is an expression for some of the components of the momentum in terms of the 
variables of the reduced phase space) 

rAoMo) = dr l-VBP^Aoir,0o) + lcZ'\r,do)^ (61) 

and 

1 r ’"o 

fr{ro,0o) = 2 dr{{p^^qAB,r){r,9o) + G^’^*'\r,0o)) 

+ drVA{q^^{r,0o) £ dflvcP^Bif,0o) - ic'r“(A0o)j). (62) 

The Hamiltonian governing the dynamics of the reduced phase space is again given by equations (I58|> 
and ( 1551 ) : however, the expressions for ffa should now be taken from the above two formulas. This 
Hamiltonian corresponds to fixing the shift vector field to the one presented in Appendix lA 41 
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D. Where should we measure from in the full theory? 

In this section, we will compare the two implementations of the radial gauge without assuming spherical 
symmetry, discussing their problems and advantages. 

Taking into account the asymptotic behavior of the fields discussed in Appendix IB 21 we see that the 
derivation presented in Section fiV B I runs into problems. In particular, (1531) and (I56p are not well-defined 
because the integrands are, in general, not integrable at infinity. Those expressions can be viewed as 
solutions of the vector constraint 

drp^^ = -Vbp^a + (63a) 

drP^r = \qAB,rP^^ - dAp\ + (63b) 

Applying the asymptotics from Appendix IB 21 we see that DbP^ a general, be finite for r going to 

infinity and hence, it is not integrable in this limit. There are at least two possible ways to deal with this 
problem. Firstly, one can try to impose a more stringent fall off behavior of the momentum field. This is 
problematic because of the necessity of such an imposition to be dynamically consistent, i.e., preserved 
by the Hamiltonian. Secondly, one can try to replace the equations for p^a with equations for p’'“. It 
leads to the following form 

o rA , AB rC BA , ^ AB(n* \ 

OtP +g qBC,rP =-VbP +-q Cb , (64a) 

drP^'- = ^qAB,rP^^ " OaP^^ + • (64b) 

The right-hand side of the first equation is integrable in r at infinity; however, the differential operator 
on the left-hand side is much more complicated. For practical purposes, this operator would need to be 
explicitly inverted, so that the dependence of the solution on the reduced phase space variables is known. 
These problems are the drawbacks of the “measuring from zero” implementation of the radial gauge. For 
some applications of the presented construction, it may be an advantage of this implementation that, as 
shown in Appendix I A 31 it can be guaranteed that the shift vector Nh belongs to the generators of Diffobs- 
Moreover, the problems at infinity may be cured by fixing the canonical data at some finite boundary 
and implementing the construction in the bounded region only. 

The implementation of the radial gauge coined “measuring from infinity” is favorable for a few reasons. 
Firstly, unlike the other implementation, this one is in agreement with the asymptotic flatness require¬ 
ments as spelled out in Appendix [BJ Secondly, it is the one which allows for a natural description of a 
Schwarzschild black hole (see Appendix [C|) . Note also, that this solution was the one already suggested 
in i. 


E. Summary of the case without spherical symmetry 

In this section, we have derived two reduced phase spaces for general relativity. Each of them is 
parametrized by a one parameter (r is the parameter) family of intrinsic geometries of the surfaces of 
constant radial distance from the center, described by qAB, and their momenta p^^ (possibly also by 
matter fields, in case they are present in the theory). Their evolution is generated by the Hamiltonian 
constraint (155)) (where ffa are given by (1551) and ((551) or in the other case, by (15T1) and ((551) 1. Details 
concerning imposing asymptotical flatness in those reduced phase spaces are described in Appendix IB 21 


V. CONCLUSION AND OUTLOOK 

In this paper, we have used the radial gauge to construct reduced phase spaces for general relativity 
with and without assuming spherical symmetry. 

The spherically symmetric setting allows for a straightforward treatment, since the tensorial structure 
is considerably simpler (e.g., intrinsic geometry is described by just two functions, the shift vector is a 
single function). A drawback of the construction is that it results in a nonlocal Hamiltonian, leading to 
nonlocal equations of motion. On the other hand, the construction for the spherically symmetric case 
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exposes various nontrivial features of the idea and therefore may serve as a toy model for the full case. 
More importantly, it suits very well a certain quantization scheme (discussed in @ and M) , where it 
plays an important role in the definition(s) of the spherically symmetric sector of the theory. This is 
of particular interest, since up to date, the treatments of the simplest (spherically symmetric) collapse 
scenarios, including quantum gravity effects, have only been performed in midisuperspace models. 

The reduction of the nonsymmetric case leads to the same kind of nonlocalities in the Hamiltonian. 
However, the reduced data still possess a clear geometrical interpretation and allow for a quantization of 
the system. Furthermore, the splitting of the description into surfaces of constant radial distance from 
the center and the radial direction is reflected in the quantum theory in a remarkable simplification. The 
geometry is effectively described by a one parameter family of 2 + 1 dimensional geometries, opening the 
possibility of formulating computable models of quantum dynamics. 

There are a few ways in which the presented formalism can be developed further or applied: 

1. The reduced phase spaces suit well a quantization using the loop quantum gravity techniques. 
Analysis of the symmetric and nonsymmetric cases leads to a definition of a reduction to spherical 
symmetry on the quantum level. In fact, different definitions of spherical symmetry can be given, 
retaining different sets of degrees of freedom of the full theory. For the details, we refer the reader 
to 0 and [l^ . 

2. An advantage of the spherically symmetric setting is that radial geodesics spanned from the central 
point never cross each other. Generically, this does not happen if space is not symmetric. Therefore, 
in general, the construction we introduced in the nonsymmetric case will break down due to the 
formation of caustics of the radial geodesics. We excluded this by demanding that the coordinate 
system we use is global. In general, this issue may also be addressed by modifying our construction, 
so that we require the radial gauge condition to hold only in some neighborhood of the central point. 
To be able to consistently perform this modification one would need a (delicate) introduction of a 
boundary up to which the gauge condition holds. 

3. An interesting question to ask is whether the current construction can be generalized to incorporate 
spatially compact settings. It seems that at least the spherically symmetric case is open for such 
a generalization, since the point antipodal to the center of symmetry will again be a center of 
spherical symmetry. 

We leave the development of the ideas presented in the second and third points for future research. 


Appendix A: Shift vector fields corresponding to the chosen gauges 

Our aim in this appendix is to compute the specific shift vector fields, we call them Nh, which, when 
included in the Hamiltonian, correspond to the radial gauge in the variants presented above. 


1. Spherically symmetric case — measuring distance from zero 

The Hamiltonian given by m is the Hamiltonian of the reduced phase space obtained after fixing the 
gauge in the variant in which we “measure from zero”. The question we address in this section is: What is 
the shift vector field defined in the context of the unreduced phase space such that the Hamiltonian from 
that phase space (including the vector constraint smeared with the shift we seek) preserves the gauge we 
chose? In other words, we want to find Nh such that 

H[N]\,,u^,_i,^ = H[N] + C[NH]. (AI) 

This can be easily done by substituting m instead of (II5cp into the Hamiltonian constraint and sepa¬ 
rating the appropriate terms. After a few manipulations, we can read off 



Note, that this shift vanishes at zero as it should for regular spacetimes in the context of spherical 
symmetry. 
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2. Spherically symmetric case — measuring distance from infinity 

The shift vector which corresponds to the Hamiltonian given by (1261) can be obtained in the same 
manner yielding 


Note that the above shift is not automatically vanishing at zero. This fact will be exploited in Appendix 

o 


3. Full theory case — measuring from zero 


In this section, we will discuss the shift vector field which corresponds to fixing the radial gauge in the 
context of the full theory. The strategy is the same as we used in the context of spherical symmetry. As 
a starting point, we use 

p'’airo,6o) = P'^a{ro,0o) - C[ M ], (A4) 

a{ro:do) 

where 

r M ir,9) = eir-ro)Sie,eo)dAo, 

I Ao(ro,6'o) / A r\ 

I M (r,9) = ^e{r-ro)S{e,eo)dr-^Q{r-ro) rdfq^^(f,9)(dAS{0,eo))dB, 

Kr{ro,9o) 

in the definition (1581) of the Hamiltonian. Then separating the Hamiltonian constraint part and extracting 


the vector constraint part, we can read off the relevant shift vector field Nh, such that 

H[N]\,,u^,_i,^ = H[N] + C[NH]. (A6) 

It turns out that 

N^Hir,e) fjdro + qABP^^ + M) (A7a) 

N^{r,e)=£ dro [q^^ (-Vs + 2^^)) (ro, d) (A7b) 

r^oZ/V l\ 

droq^^{ro,0)dB df l-j=={-p^’'+ qcoP^^ +-U)] (A7c) 

where we used the shorthand notation 

ta{ro,eo)--=C[ M I (A8) 

a{ro,eo) 


A thing to notice is that the above field is almost of the kind spelled out in ([55|) if one identifies 

COrr = -£= (-(?'■’■ + qABP^^ + ), (A9a) 

Vdet q 2 

N 

LOrA = —=={-^p'" A + 2tA)- (A9b) 

Vdetg 

As it was pointed out in Q, the radial gauge we are using does not fix the spatial diffeomorphism gauge 
completely. In fact, it fixes diffeomorphisms called Diffobs, leaving a finite (six) dimensional, residual 
diffeomorphism freedom. Therefore, it may be desirable to stay within the Diffobs class when reducing the 
phase space. The part missing from (IA7I) is the contribution from uj at zero (compare with equation (1361) 1. 
However, it turns out that we can consistently fix that problem. Notice, that the required contribution 
from w at zero vanishes when we are in a point of the constraint surface in which 
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Let us postpone the discussion of the meaning of that condition for now, to notice that the Poisson 
bracket of it with the gauge fixed Hamiltonian is again proportional to such a term. This means that if 
we require our initial data to satisfy that condition, the evolution we have formulated will preserve it. To 
understand the geometrical meaning of the condition we just imposed, let us rewrite it in terms of the 
extrinsic curvature 

lim|^-^(-p- + q^BP^^)j = {NK„) (0). (All) 

Requiring that it vanishes is equivalent to requiring 

{NKjj) (0) = 0 V7,j, (A12) 

but that means 

-igrj(0) + 5(,iVj^(0) = 0. (A13) 

Since we know that in the coordinates I we have qij{0) = 6ij constantly in time, we obtain a condition 

%iVj^(0)=0. (AM) 

This condition can easily be interpreted in the light of Q- It just means that from the initial data 
allowed by the construction without the condition (jAlOII , we chose such that no relative deformations of 
the directions associated with the central observer take place as the system evolves, or equivalently, such 
that the central point (and its infinitesimal neighborhood) is flatly embedded in the spacetime generated 
by the evolution. 

The remaining task is to express the condition (lAlOl) in the reduced phase space variables. In the 
reduced phase space, we have 


1 

/^(ro,0o)=-- / dr{{p^^qAB,r){r,eo) + Cr^\r,0o)) 

^ Jro 

X ooQ/ 1 \\ 

^ I df\VcP^B{f,9o)--CT^\f,9o)\\, (A15) 


so we can rewrite the limit 
N 


lim, _ 

r^o ^det q 


(-/'■+p^^gAB)l = 


) 


+ 5/." + cr")-jr dfV^ (,-»»(f) J,- (VcP^B - 


Both the nominator and the denominator vanish as r ->• 0, so using the L’Hopital’s rule we get 


lim 


N 


^_prr^pAB^AB)] = 


0 \ Vdet q 


= 7V(0)lim 

r-»-0 


{p^^qAB) 


1 ^AB ^ 

2P qAB,r 


+ VA{q^^ C{VcP^ 


1 /^matt 
2^B 


)) 


(Vdetg) 


(A17) 


Now we need to look at each of the terms separately 

p'^^qAB-r^ so {p^^qAB) r, 

P"^^qAB,r ~ r, 

C matt 2 

r " ^ ^ 

■DcP^B-lcT^^-r^ so VA^q^^ £^lvcP^B-lcT']]-r^, 
\/det q^ so (Vdetg) ^ ~ r. 


(A18a) 

(A18b) 

(Al8c) 

(A18d) 

(Al8e) 
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therefore, using the L’Hopital’s rule again, we get 


lim I 


N , rr , AB ^ \ HH 


\/detq 


i-p^^+p^^qAB)] = N{0)lim 


{p^^qAB) - ^{p^^qAB,r) 


r-^0 


(Vdetg) 


(A19) 


and the right-hand side is an expression involving only the reduced phase space variables. 


4. Full theory case — measuring from infinity 

In this section, we present the shift vector Nh which corresponds to the Hamiltonian from Section 


IIV Cl Rewriting the expressions from (ICT) and (1621) in the form 

ra=p\-C[Ml (A20) 

a 

we can express the reduced Hamiltonian as 

H = H+C[Nh], (A21) 

with 

droi-^=={-p^'' + qABP^^ + -tr)Uro,d) (A22a) 

N^ir,0) =-dr 0 (^q^^ b + 2 ^ 5 ) j (ro, 0) (A22b) 

1 Z'00 ^00 / J\f 1 \ 

-2 Jr droq'^^{ro,d)dB dfi-^=={-p^'^ + qcDP^^ + -tr)Ur,0), (A22c) 


where ta are terms proportional to the vector constraint. This shift vector vanishes at infinity; however, 
its behavior at zero is more obscure. 


Appendix B: Asymptotic conditions on canonical data 

In this appendix we will discuss the falloffs of the canonical fields for which the treatment presented in 
the main part of the paper is well-defined. We will restrict the discussion to the gravitational sector of the 
phase space, since it is the nontrivial one. Our analysis will be based on the treatment of asymptotical 
flatness presented in [Hi- We demand that there exist asymptotically Minkowskian coordinates (a;*) 
(additionally, we define a coordinate p := \/Sij and also, n* := in which the fields behave as 

q^jix'^) = (HI) 

P 

= + (B2) 

P 

where, as the notation suggests, Sij and are functions of the angles only, and they satisfy the so-called 
parity properties 

Sij{n'") = Sij{-n'"), (B3) 

= (B4) 

while by we denote terms which vanish at infinity as their first derivatives vanish as 

^-(”■+1)^ and so on. The asymptotical behavior of the lapse and shift is then 

Nix^) = k{n'^) + 0°°ip-^), 

N\x^) = k\n^) + 0°°ip-^), 


(B5) 

(B6) 
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where 

kin^) = -k{-n^), (B7) 

k\n^) = -k\-n^). (B8) 

For such fields, one can check that all the integrals in the canonical theory are convergent “at infinity" and 
functionally differentiable. As was argued in (l^ . dropping the parity properties of the lapse and shift 
gives rise to the ADM four momentum. Here, we loosen the parity property of the lapse only; therefore, 
the asymptotics of lapse functions we consider are 

Nix'^) = l + k{n^) + 0^{p-^). (B9) 


1. Spherically symmetric case 

The translation of the asymptotic conditions spelled out in the preceding paragraph to the variables 
used in spherical symmetry was carried out in |12l |. It gave 


A(p) = l + - + 0“(p-i-), 
p 

(BlOa) 

i?(p) = p+0'-(p-^). 

(BlOb) 

Pa(p) = 0“(p-), 

(BlOc) 

Pfl(p) = o“(p-i-). 

(BlOd) 

iV(p) = l + 0“(p-^), 

(BlOe) 

N^ip) = 0“(p-), 

(BlOf) 


where it is easy to check that the formula (|5]) for the mass functional, yields m in the limit of large p 
when the above asymptotics are employed. This justifies the usage of the same symbol in the definition 
of the functional and in the asymptotic considerations. It should be noted that the shift may be allowed 
to be finite at infinity, since it corresponds to an asymptotically odd shift in agreement with (IB8I) . 

Such conditions require a boundary term to be added to the Hamiltonian constraint to secure its 
functional differentiability. That boundary term is the ADM energy that in this context is given by 

E = lim -p{lS? + ^ -= lini —p(A^ - 1) = m. (Bll) 

P-.00 2 p^ p p-foo 2 

Note that since the beginning of the discussion of asymptotical flatness we used the variable p to denote 
the radial coordinate denoted by r in the original papers (Em). This is because in the current paper, 
the label r is reserved for a specific radial coordinate, namely, the one that measures the proper spatial 
distance. Since in our analysis we want to use the variable r, we need to transform the conditions (IBIOI) . 
What is the relation between r and pi Making use of the fact that A is a density and knowing on one 
hand that it satisfies (IBlOal) and on the other that it is equal to 1 as a function of r, we find the condition 

A{p)dp=l-dr. (BI2) 

After integrating this condition and finding the leading terms in the expansion of p as a function of r, we 
find 


p{r) = r - mlogr + c+ 0°°{r 


(BIS) 


where c is a constant whose role is to fix the specific value of p at a given r. For simplicity, we choose 
it to be zero in the following. Using this relation in (IBlOll yields the asymptotic conditions on the fields, 


which we use in the current paper, 

A(r) = 1, (B14a) 

R{r) = r - mlogr + 0°° (B14b) 

PA(r) = 0“°(r-^), (B14c) 

= (B14d) 

Af(r) = l + 0°“(r-"), (B14e) 

iV’'(r) = 0°“(r-'^). (B14f) 
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Note, that condition (IB14a|) is satisfied not only in the limit of large r, but for all its values. 

Having defined the asymptotic behavior of the fields, we should address the problem of boundary 
terms. It turns out that the vacuum Hamiltonian and vector constraints are both well-defined (meaning 
integrable at infinity). However, to ensure the functional differentiability of the Hamiltonian constraint 
we have to add a boundary term analogous to (IBllI) with 

1 d2 O D JDf 

E = lim -r(l + ^- = m. (B15) 


2. Full theory case 

In this section, we will specify the asymptotic behavior of the fields used in Section ITVl It amounts to 
a translation of the conditions we are using into the variables we work with. 

We start by rephrasing the condition (IBII) in spherical coordinates such that cc® = pn^(0) (note that 
those coordinates differ from the ones we want to use eventually by a rescaling of the radial coordinate). 
It reads, in particular, 


Qpp = 1 + —n^n^Sij + 0°°{p ^ '^), (B16a) 

P 

QAB = P^VAB + pn\n]gSij + 0°°(p^“''), (B16b) 

where pab is the metric of a unit sphere. To switch to the variables we want to use, we need to find r 
such that qrr = 1. In order to do that, we compare lengths of intervals of radial geodesics and obtain the 
relation 

/ = (B17) 


which leads to 


r = p+^rPn^Sij\ogp + c + 0°°{p '^). (B18) 

Like in the previous section, the role of the c is just to fix the relative values of r and p, so we choose it 
to be zero in what follows (it can be restored easily). Moreover, since the angular variables coincide, we 
can use instead of rP from now on. Inversion of (|B18I) gives 

p = r - rv' sijlogr + 0°° (B19) 

The asymptotics of the metric can now be given explicitly 

Qra — ^rat (B20a) 

QAB = -n^n'^Sijrlogr)pAB +rn^^A'r>''!BSiJ + (B20b) 

Transforming the momentum is more involved, because it is a tensor density. Using the notation 
^(0) ■= ; 3 det(§^) and ff{0) := we find 

p\ = n{0)ni{0)n-’{0)t^j{0) + 0°“(r-'^), (B21a) 

/a = rn{0)ni{0)n-^Ai&)t'Ji&) + (B21b) 

p^B = n{0)ff{0)n-^BWt' jW + 0“(r-^)- (B21c) 

The fact that the leading terms in the last and second to last lines are not vanishing for large r is the 
source of the problems with the “measuring from zero” implementation of the radial gauge discussed in 
Section IIVDI In the “measuring from infinity” implementation, however, such conditions render all the 
integrals well-defined. 
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Having the asymptotics of the metric at hand, we can compute the ADM energy. We obtain 

E = lim ,(/ ^K^Jdei{qAB)(^^0 - J ^K^det(qAB)d^0) (B22a) 

= -^ y' (n^rf^ - sjj det{q)(f 9, (B22b) 

where in the first line is (the trace of) the extrinsic curvature of the surfaces of constant r as embedded 
in the spatial slice E, while the ° symbol denotes a flat contribution which needs to be subtracted for the 
final result to be finite [2^ [2l| (see also for a discussion in similar coordinates). 


Appendix C: Generalization of the radial gange to spherically symmetric spacetimes singnlar at 

zero 

Although the “measuring from infinity” variant of the reduced phase space construction with the asymp¬ 
totical behavior of the fields described above seems well suited for a treatment of Schwarzschild black 
holes, in fact, it is somewhat problematic. The Schwarzschild solution is expressed in variables of [T^ 
by“ 

A(p) = — R{p) = p, Pa(p) = 0 = Pr(p). (Cl) 

P 

To find expressions for the canonical fields in the radial gauge, we use again the condition (IB12II . Inte¬ 
grating it (for concreteness we choose r to measure the proper distance from the horizon), we find 

+ (C2) 

V p \Tn my p / 

To identify the Schwarzschild solution we should now invert this function, finding p{r). Then 

A(r) = l, R{r)=p{r), (C3) 

would be the Schwarzschild solution in the radial gauge. Unfortunately, inverting that function explicitly 
is very hard.^^ It can be checked that the inverse function fulfills the asymptotic behavior spelled out in 
(jB14bp . but using it to describe the Schwarzschild solution is highly impractical. Therefore, it is desirable 
to modify the construction slightly in order to accommodate the Schwarzschild solution more easily. 

Being guided by that aim let us recall that Schwarzschild metric expressed in Gullstrand-Painleve 
coordinates has the form 


ds"^ =-{I-—)dt^+ 2x[^dtdr + dP + Pdn^. (C4) 

r \ r 

The idea of using those coordinates in the context of the canonical analysis of spherically symmetric 
general relativity is not new (see, e.g., (23 - [^ and references therein). They are particularly useful for 
us, since we see that on t = const surfaces the variable r measures the proper (spatial) distance. Note 
however, that the shift vector read off from the above form of the metric is 


AT(r) = 



(C5) 


which means it does not vanish as r goes to zero, contrary to the case for regular spacetimes considered 
in Section uni It does vanish, however, when r goes to infinity. It is this behavior of the shift that we 


Note that in equation (58) of 0, only the leading orders of the Schwarzschild solution are given. Hence, they satisfy 
the Hamiltonian constraint only in the limit of large p (or Irl in the notation of 0)- 

The problem can be reduced to inverting the function x = (fyl + for y{x), which is a modified defining problem 

of the Lambert W function (see 0 ). 
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will incorporate in the present modification of the radial gauge construction. Because of this behavior 
we will work in a setting similar to the “measuring from infinity” variant presented in Section [ill PI 
From the form of the metric presented in (IC4I) , we can see that for the Schwarzschild spacetime 


A(r) = 1, A(r) = 0, (C6) 

R{r) = r, R{r) = 0, (C7) 

N{r) = 1, 7V”(r) = (C8) 

Using the definitions of canonical momenta spelled out in ([4]) (with the convention y = 1), we obtain 

Pnir) = (C9) 


Inspecting equation (1261) , we see that the reduced Hamiltonian constraint vanishes for the above data. 
The shift vector preserving the gauge expressed in (j27ll is equal to the one read off from the Schwarzschild 
metric, as written down in (IC5I) . Finally, computing the equations of motion (with the x = 1 convention) 


—i?(r) 
N ^ ^ 

—Pni-r) 

N ^ ^ 



Fjr) \ 

i?2(f) ) 


-R"{r) + 


PUr) 

R{r) 


PR{r)F{r) 

R^{r) 


F\r) 

R^{r) 



F{f) \ 
R^{f)) ’ 


(ClOa) 

(ClOb) 


one can verify they are satisfied, leading to the conclusion that a Schwarzschild spacetime is indeed a 
solution of the equations of motion of the reduced phase space presented in this appendix. 

An apparent problem, which arises in this treatment, is that the fields in equations (IC6|1 to (|C9ll do 
not satisfy the falloff conditions spelled out in Appendix [Bj This change of asymptotics of the fields can 
be understood as stemming from a different choice of slicing in the corresponding spacetime. The falloff 
behavior of the fields in the current context can be obtained by inspecting what sort of variations of the 
canonical variables are allowed in order to maintain well-definiteness and functional differentiability of the 
constraints. It yields the asymptotic behavior which needs to be employed in the phase space presented 
in this appendix (and propagating in a clear-cut way to the asymtpotic behavior of the corresponding 
reduced phase space), namely. 


A(r) = 1, 

(Clla) 

R{r) = r + 0^{r-^), 

(Cllb) 

PA{r) = ^2mr + 

(Cllc) 


(Clld) 

N{r) = \ + 0°°{r-^), 

(Clle) 

Ar’'(r) = y^+0~(r-5-^). 

(Cllf) 


Additionally, a boundary term equal to the Schwarzschild mass needs to be added to the vector constraint 
to ensure its functional differentiability, replacing a boundary term for the Hamiltonian constraint present 
in the previous asymptotic analysis. 
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